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ABSTRACT 

MUSE observations of NGC5813 reveal a complex structure in the velocity disper¬ 
sion map, previously hinted by SAURON observations. The structure is reminiscent of 
velocity dispersion maps of galaxies comprising two counter-rotating discs, and may 
explain the existence of the kinematically distinct core (KDC). Further evidence for 
two counter-rotating components comes from the analysis of the higher moments of 
the stellar line-of-sight velocity distributions and fitting MUSE spectra with two sep¬ 
arate Gaussian line-of-sight velocity distributions. The emission-line kinematics show 
evidence of being linked to the present cooling flows and the buoyant cavities seen in X- 
rays. We detect ionised gas in a nuclear disc-like structure, oriented like the KDC, which 
is, however, not directly related to the KDC. We build an axisymmetric Schwarzschild 
dynamical model, which shows that the MUSE kinematics can be reproduced well with 
two counter-rotating orbit families, characterised by relatively low angular momentum 
components, but clearly separated in integral phase space and with radially varying 
contributions. The model indicates that the counter-rotating components in NGC 5813 
are not thin discs, but dynamically hot structures. Our findings give further evidence 
that KDCs in massive galaxies should not necessarily be considered as structurally or 
dynamically decoupled regions, but as the outcomes of the mixing of different orbital 
families, where the balance in the distribution of mass of the orbital families is cru¬ 
cial. We discuss the formation of the KDC in NGC 5813 within the framework of gas 
accretion, binary mergers and formation of turbulent thick discs from cold streams at 
high redshift. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: formation - galaxies: 
evolution - galaxies: kinematics and dynamics ~ galaxies: nuclei - galaxies: structure 


1 INTRODUCTION 

The vast majority of early-type galaxies are kinemati¬ 
cally simple objects. Their stellar motions can easily be 
predicted from their light distributions, confirming they 
are largely oblate anisotropic systems llCappellari et all 

* E-mail:dkrajnovic@aip.de 


| 2013b|). The SAU RON (Ide Zeeuw et al.l|2002ll and ATLAS®® 

( Cappellari et all |2011^ surveys showed that within their 
half-light radius the majority of early-type galaxies (~ 85 per 
cent ) have simple disc-like velocity maps a ligned with their 
light distributions llKrainovic et al.l l2008l . 1 201 il l , and that 
they have relatively high specihc angular momentum, and, 
there fore, are referred to as fast rotators llEmsellem et al.l 

iml). 
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The topic of this paper is about an exception to the rule. 
The minority of early-type galaxies , those with irregular ve¬ 
locity maps llKrainovic et al.ll201ll i and lo w specific angu¬ 
lar momen ta, referred to as slow rotators jEmsellem et al.l 
l2007l.l201ll\ can also exhibit spectacular kin ematic features 
llEmsellem et aL|l2004l : iKrainovic et al■ll2011^ . The most eye¬ 
catching features on velocity maps are certainly kinemati¬ 
cally distinct core^ (KDC), which come in a variety of forms. 
Often they are found in objects which have no net rotation 
at larger radii and only the KDC part of the galaxy rotates. 
In other cases, when the outer body shows rotation, the axis 
of rotation of the KDC and the outer body are misaligned 
by approximately 90°or 180°(counter-rot ation). Other mis¬ 
alignments are also ob served, but are rare dFranx et al.lll99ll : 
IKrainovic et al.ll201ll ). 

The fi rst object with a KDC to be discovered was 
NGC5813 (lEfstathiou et al1ll980l . [19821 ). the subject of this 
study. As kinematic studies of elliptical galaxies became 
more numerous in the late 1980s and early 1990s, the number 
of known KDCs also increased llFranx fc Illirigw orthl Il988l: 
Jedrzejewski_&_^checht^[W8^ iBendeil 19881 : iFranx et ^ 


19891 : Carollo &; DanzigCT 1994) . From the earlier studies 


it was clear t hat KDCs are frequent in massive ellipticals 
llBendeilll988l) . but the ATLAS®° survey showed that, while 
the KDC host galaxies span a range of masses, they make 
only a small fractio n (7 per cent) of the t otal population of 
early-type galaxies (IKrainovic et al.ll201ll) . KDCs are, how¬ 
ever, numerous among slow rotators (~ 4 2 per cent), but 
very rare among fast rotators (~ 2 per cent: lEmsellem et al.l 

M). 

Detection of KDCs is obv i ously a question of spa¬ 
tial resolution. iMcDermid et al] ll2006l ). using OASIS, an 
integral-field spectrograph with higher spatial resolution 
than SAURON, found a number of small KDCs, typically of 
only a few hundred parsecs in diameter. These KDCs were 
actually found in fast rotators and seem to be only counter¬ 
rotating (oriented at 180° with respect to the main body). 
The same study showed that these small KDCs should be 
considered different from the large KDCs found (mostly) in 
slow rotators, as they have younger stellar ages, often asso¬ 
ciated with gas and dust. Moreover, the authors argued that 
the small KDCs are essentially the product of recent accre¬ 
tion of relatively minor amounts of gas, which rejuvenates 
the nuclei, and the KDC is visible temporarily only because 
the mean velocity is a luminosity weighted measurement, 
and as such is biased by bright young stars. As the stel¬ 
lar population of the small KDC stars ages, its imprint on 
the observed kinematics on the host body will diminish and 
eventually disappear. 

The large KDCs are made of old stars, with colour 
or stellar population gradients similar to other early- 
type galaxies of the s ame mass dFran x fc IllingworthI 19881: 
Caro]lo_e£^jJjl997a|^ Davies et al.l 1200ll : iKuntschner et al l 


20101 : McDermid et al. 20151 ). but their origin is not 


pletely understood. The distinct nature of the angular mo¬ 
mentum vector of the KDC with respect to the rest of 
the host galaxy suggests a KDC of external origin, and 


^ Sometimes they are called kinematically peculiar cores, or kine¬ 
matically decoupled components or cores, but the decoupling na¬ 
ture of these structures is not confirmed. 


an evidence for mergers. The observation that the sur¬ 
face brightness pr ofile of NGC 5813 can be described with 
two componen t s (lEfstathiou et ^ Il980l ). was generalised 
by i ormendtl lll984h . who suggested that the KDC in 
NGC 5813 is a remnant of a small galaxy, and a similar 
argument could be applied to other KDCs. The light pro¬ 
files of KDC host galaxies are, however ^ not different from 
other slow-rotators of similar masses dPorbes et al.l 1 19951 : 
ICarollo et al.l[l997al f9). and the similarity of stellar popu¬ 
lations within the KDC and the outer body of the galaxy 
is not consistent with t his ’core-within-a-core’ scenario 
(iFranx fc Illingworthlll98^ . In particular, if cores originated 
as small ellipticals, one should observe positive metallicity 
gradients (core would be metal poorer than the rest of the 
gala xy), but the opposite to thi s is observed in KDC hosts 
(e.g. Efstathiou fc Goreaj 198^ FYanx&IlhngworD^ 198^ 


Bender fc Surmal Il992l: iMehlert et al.l Il998l: iDavies et all 


200ll : Kuntschner et ani201(]| ). 


Other possible scenarios for the formation of KDCs in¬ 
clude the disruption of a satellite galaxy during a minor 
merger, where the angular momentum of the KDC is the 
consequence of the orbital a ngular momentum of the satellite 
(IFranx fc Illingworthlll98§ ). a dissipative merger of gas rich 
(spiral) galaxies, or accretion of external gas, where in both 
cases some of the acquired gas settl es into the nucleus of th e 
main progenitor and forms stars (iBender fc Surmi3 Il992l) . 
While KDCs seem to be indica tors of past merger history, 
Ijedrzeiewski fc Schecht^ (1 19881 ) pointed out that photomet¬ 
ric evidence (i.e. shells) are not often associated with KDCs. 
This, however, is consistent with an early formation (sup¬ 
ported by the old ages of stars within KDCs), where galaxies 
are built through hierarchical merging of small bodies, pos¬ 
sibly partially gaseous, at an early epoch dBender &: Surmal 
1 19921) . 

The formation of KDCs has to be considered within the 
cosmological picture of galaxy formation, namely that there 
seems t o be two phases o f mass assembly for massive galax¬ 
ies (e.g. lOser et aLll2010|) . where the first one i s dominated 
by vigorous star formati on fed by cold streams dKeres et al.l 
I2OO5I : IPekel et al.l l2009l ). followed by the accretion of al- 
ready formed stars f r om satellites (e.g. iNaab et al.l l2007l: 


„ __ _ !. g. IlNa; 

Feldmann_et_aL I I2OIII : Ijohansson et al.ll2012l : 

2 OI 2 I) . 


Lackner et al.l 


There is one other group of kinematically complex 
galaxies, which are found on both sides of the slow -- fast 
rotator division line, but make up only 4 per cent of the to¬ 
tal population of early- types in the nea r by U niverse. These 
galaxies were named bv iKrainovic et aP d201ll ) double-sigma 
galaxies, or 2a, because of unusual velocity dispersion maps. 
Velocity dispersion maps of 2a galaxies are characterised by 
two symmetrically off-centered peaks along the major axis. 
The velocity maps can, however, have a variety of structures: 
from no net rotation to counter-rotating cores, and even a 
multiple reversal of rotation direction (for examples of veloc¬ 
ity and vel ocity dispersion maps of 2a galaxies see figs G.5 
and C.6 of IKrainovic et ^l201ll ). 

The 2a galaxies are simpler systems than KDG hosts, 
and their internal structure is better unde rstood. They are 
also ty pically less massive (< 5 x Mr;^. ICappellari et al.l 
l2013bl ) than galaxies with KDCs, such as NGC 5813, but 
KDCs are also found in galaxies populating the same range 
in mass as 2cr galaxies. The two-dimensional spectroscopy 
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of these galcixies revealed a double-peaked structure in the 
line-of-sight velocity distribution (LOSVD) indicating they 
are made of two separate ki nematics components which 
counter-rotate (e.g. NGC 4550, iRubin et al.lll99^ : lRix et alj 
I 1992 I I. Dynamical models based on the orbit superposition 
have shown that such galaxies are indeed made of two 
counte r-rotating components, both having high angular mo¬ 
menta (|CaDDellari et al1l2007l l. Furthermore, the two compo- 
nents are characterised by different stellar popu lations (e.g. 
ICoccato et al.|[201ll . I2OI3I : I Johnston et al.|[2013l l. Given the 
flatness of these galaxies, and the high angular momentum of 
each separate component, 2a galaxies are consistent with be¬ 
ing made of two counter-rotating discs of comparable masses, 
where the details of their kinematic features (rotation or no 
rotation in the velocity maps, the shape and the separation of 
the two peaks in the velocity dispersion maps) depend on the 
relative masses and characteristic sizes of the discs. For ex¬ 
ample, NGC 4550 seems to be made of equally massive discs, 
while the counter-rotating discs in NGC 4 473 are approxi¬ 
mate ly split 30-70 per cent in mass (e.g. ICaPDellari et al.1 
l2007ll . 

The origin of 2cr galaxies is currently associated with 
accretion of gas, where the accreted external gas settles in 
the equatorial plane, but with an opposite angular momen¬ 
tum. The source of the gas can be a nearb y object, such 
as in the case of NGC 5719/NGC 5713 pair llVergani et al.l 
l2007l : ICoccato et al.|[2^ , or, as recent numerical simula¬ 
tions suggest, a 2a galaxy could sit in a locus of two cold 
streams, which feed the galaxy with gas and exert different 
torques resulting in the formation of the prograde and ret¬ 
rograde disks of relative contributions propo rtional to the 
properties of the streams llAlgorrv et al.ll2014l l. 

In this paper, we study NGC 5813 with MUSE (Bacon 
et al. in prep.), focusing on its KDC, and the fact reinforced 
by new observations that its KDC is the result of the super¬ 
position of two counter-rotating structures, similar to what 
is found in 2a galaxies. This has strong implications on the 
formation history of this kind of galaxies and the origin of (at 
least some) KDCs in massive slow rotators. We argue here 
that the KDC in NGC 5813 is indeed an orbital composite, in 
a simi lar way as found in NGC 4365 by Ivan den Bosch et al.l 

II2OO8II . 

We start with a summary of general properties of 
NGC 5813 (Section [2|) relevant for this study. The obser¬ 
vations and data reduction are presented in Section [31 The 
main evidence supporting our picture are the MUSE stel¬ 
lar kinematics presented in Section [4] and the Schwarzschild 
dynamical models from Section [S] In Section [5] we present 
emission-line kinematics from MUSE obeservations, which 
show that the present gas distribution is not related directly 
to the KDC. We summarise the main results and discuss 
the possible origin of the KDC in NGC 5813 in Section [T] 
Throughout this paper we assum e NGC 5813 is at a distance 
of 31.3 Mpc llTonrv et al.ll200lll . 


2 GENERAL PROPERTIES OF NGC 5813 

NGC 5813 is one of the two b r ightes t galaxies in the 
NGC 5846 group (iMahdavi et al.l I2OO5I I. This is an iso¬ 
lated group which consists of three parts dominated 
by NGC 5846, NGC 5831 and NGC 5813. Spectroscopic 


observations jEfetathlouet^ IRSd; Davies et al.l 1 19831 : 


ICarollo et al.l Il993l : lEmsellem et al.l l2004l l revealed that 
NGC 5846 shows very low mean motion, while NGC 5831 and 
NGC 5813 both harbour KDCs. NGC 5846 and NGC 5813 
have simil ar K-band magnitude (-25.01 and -25.09, re¬ 
spectiv ely, ICappellari et al.l 1201 ill , while ICappellari et ah] 
(l2013bll derived masses of 3.7 and 3.9 x IO^^Mq, respectively 
(constrained within the SAURON observations). NGC5831 
is a smaller and a less massive system. 

The optical imaging and spectroscopy show that 
NGC 5813 is unusual both kinematically and photometri¬ 
cally. The light profile in the observations of lEfstathiou et aP 
(I 1982 I I indicated a "two-component structure" and was 
markedly different from other elliptical galaxies. There are 
a number of ways one could parameterise such a light dis¬ 
tribution, and this is evident in the variety of approaches in 
the recent work. 

The n uclear region can be fit by a Nuker double 
power la w (|Lauer et al.| Il995l l. but also with a cor e -Sersi c 
(ICraham et al.l l2003l l. as done by iRichings et al.l (1201 il l. 
Outer regions (beyond 2'.'5) are well fit with a combina - 
tion of a Sersic and an expone ntial (iKra inovic et al. 2013al l. 
while fitting the HST imaging iDullo fc GrahaiJ i 2012l l and 
iRusli et al.l ll2013h agree that the best fit is achieved through 
a combination of a Sersic and a core-Sersic profiles, where 
the core-Sersic has a low Sersic index (n=2.8 or 2, respec¬ 
tively), while the outer profile is close to exponential. There¬ 
fore, the outer profile of NGC 5813 is remarkably similar to 
tho se found in disc galaxies . This outer exponential profile 
led IDuUo fc GrahanJ ll2014h to classify NGC 5813 as an SO 
galaxy. 

The isophotes sh ow a change in both thei r flattening and 
their orientation le.g. iDullo fc Graham|[20l3l . Not consider¬ 
ing the behaviour in the sub-arcsecond region, the change 
in ellipticity is prominent beyond 10", roughly where the 
KDC stops, and changes from being approximately constant 
at 0.1 (within 10") and continuously increasing to about 0.3 
in the outskirts. The position angle exhibits a slightly more 
complex behaviour. Within the central 20", it continuously 
changes with radius (within about 15°), but there is a no¬ 
ticeable difference between the KDC region (< 10") where it 
reaches a maximum (about 150°) and the outer region where 
it fl attens to abou t 135°. _ _ 

ICaroIlo et al.l (Il997bll and iTran et al.l (I2OOII I detected 
dust in the central region of NGC 5813. The dust is ob¬ 
served in two distinct morphological configurations, as fil¬ 
aments and what appears to be a nuclear disc. The fila¬ 
ments are similarly distributed like the ionised gas converg¬ 
ing on the nucleus, while the dust disc, a thin slab of dust 
about 1 " in size, cuts the nucleus in half and it is oriented 
at ~ 148°. This orientation is on ly a few degrees diffe rent 
from the orientation of the KDC (IKrainovic et al.ll201ll l. as 
well as the nuclear disc-like kinematic feature seen in the 
ionised gas, which we present in Section 16.21 The V-I HST 
colour images indicate that what appears to be a double nu- 
cleus in NGC 5813 is l ikely a consequence of dust obscuration 
(ICaroIlo et al.|[r997bl') . 

The sub-group o f NGC 5813 is rich in diffuse X-ray emis¬ 
sion. IMahdavi et al.l (l2005ll find that X-ray emission profile 
of NGC 5813 differs from the one of NGC 5846, in the sense 
that it is less steep, more luminous and more extended. This 
could be interpreted as a consequence of different forma- 
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tion histories, where NGC5846 experienced more interac- 
tio n with smaller galaxies (for ev idence of a recent merger 
see iGoudfrooii fc TrinchierilllQQSf i . 

The AGN activity of NGG5813 is currently not high. 
It is detected in NV SS dCondon et al.lll998ti at 1.41 GHz, 
as well as at 8 GHz (iKrainovic fc Jaffell2002h . At the spa¬ 
tial resolution of the latter study (2-3”) the radio source is 
unresolved. Perhaps the most telling evidence of past, but 
also quite r ecent AGN act i vity a re the X-ray observations 
analysed by iRandall et al.l ll201ll ) , which present a case for 
three buoyant cavities, lifted by the AGN jet which operated 
in time intervals of 90, 20 and 3 million years. The cavities 
rise approximately along the minor axis of the stellar body, 
where the emission-line gas is also located, as detected on 
the Hq narrow-ban d observations ( Go^frooii et al ] 11994 
IRandall et al.l 1201 ll : I Werner et al.l 2014 ) a,n d in various 
ionised emission-line s (Section [6l Sarzi et al.l [200^ . as well 
as in lOOK cold gas ilWerner et all 2014^ . 


3 OBSERVATIONS AND DATA REDUCTION 

MUSifl is a general purpose second generation instrument 
operating at UT4 of ESO’s Very Large Telescope. It is an 
integral field spectrograph with a field-of-view of 60” x 60”, 
covered by 24 individual spectrographs, each supplied with a 
sheer of 48 mirrors arranged in 4 x 12 stacks. This configura¬ 
tion allows a wide spectral range (we use the nominal 4800 - 
9300 A), a spatial sampling of 0”2, and a spectral resolution 
ranging from R = 1500 — 3500 across the wavelength domain. 
It clearly supersedes other optical spectrographs that were 
used to observe our target. 

NGC 5813 was observed during the first commission¬ 
ing run on February 11, 2014, in fair weather conditions, 
with recorded DIMM seeing ranging between 0”6 and Cf.'8 
during the total 1.4h of integration. An a posteriori esti¬ 
mate of seeing using the ’white light’ reconstructed image 
on a few globular clusters found in the field suggests that 
the final point-spread-function (PSF) had a full-width-at- 
half-maximum (FWHM)< 0”7. The exposure sequence com¬ 
prised four on-source 20 min integrations and two offsets to 
observe a blank sky region, each of 2 min in duration, posi¬ 
tioned about 7' in the south-west direction (along the minor 
axis of the galaxy). The observing sequence was: S-O-O-O-S- 
O, where O is an object (on-target) and S is a sky (off-target) 
observation. Each on-target observation was dithered, but 
not rotated. The data were reduced using the vl.O version 
of the pipeline. 

The data reduction follows the standard steps adapted 
for MUSE and described in Weilbacher et al. (in prepj^- 
The pipeline works on pixel tables, large tables with infor¬ 
mation on all pixel values, corrected for bias and flat-field 
and their location on the CCD. Master bias, flat field and 
arc calibrations solutions are created using a set of standard 
calibrations observed on the same night. Each exposure (on- 
and off-target) was reduced separately, but in the same way. 


^ Some illustrati on of MUSE 
llBacon et al. I [20li) . 

^ A r ecent and short description 
ll2014h 


capabilities are found in 
is given in IWeilbacher et abl 


consisting of a subtraction of bias and flat-fielding. After de¬ 
termining the trace of the spectra, geometric and wavelength 
calibrations are used to transform the location of the spectra 
on the CCD into the focal plane spatial coordinates (in arc- 
sec) and wavelengths (in Angstrom). This was followed by 
the application of an astrometric solution. The standard star 
GD 108 was observed at the beginning of the same night and 
prepared within the MUSE pipeline. Notably, the response 
curve was smoothed using a piecewise cubic polynomial func¬ 
tion to remove high-frequency fluctuations. 

Sky was removed for each on-target cube separately, us¬ 
ing one of the two closest in time sky cubes. The method 
implemented in the MUSE pipeline consists of two steps. 
Firstly, the reduction software determines a representative 
sky spectrum (continuum and sky lines) from each of the 
sky exposures. These sky spectra are then adapted using the 
information on the line-spread function for each of the 24 
MUSE channels and removed from the on-target observa¬ 
tions. In particular, the first sky was applied to the first on- 
target observations, the last sky on the last two on-target ob¬ 
servations, while the middle on-target observation was sky- 
subtracted using an average sky-continuum from the two sky 
observations. 

The last step in the reduction was the merging of in¬ 
dividual observations. As there were only simple shifts be¬ 
tween the observations, the shifts were reversed and pixel 
tables were combined, producing a science-ready data cube. 
Throughout the reduction the pipeline calculates and prop¬ 
agates the formal noise at each step. 


4 STELLAR KINEMATICS 

4.1 Binning of the MUSE datacube 

Prior to extraction of the stellar kinematics, the MUSE cube 
of NGC 5813 was s patially binned using the Voronoi tes¬ 
sellation method of ICappellari fc Copinl ( 2003l L The target 
signal-to-noise ratio (S/N) was 130, where this ratio was es¬ 
timated as the ratio of the flux and the square root of the 
pipeline produced variance at 5500 A. We also created a 
binned cube with a target S/N=70 for testing and for the 
extraction of emission-line kinematics (see Section |^ . 

For the reliable extraction of stellar kinematics it is of¬ 
ten sufficient to have a lower S/N (« 40 per A), but we chose 
a relatively high S/N in order to improve on the systemat- 
ics possibly introduced by the specific dithering pattern used 
during the observations. The median S/N per pixel at 5500A 
of the MUSE data cube is 10, with the central regions reach¬ 
ing 100. The final choice of the target S/N was a compromise 
between the fact that we wanted to limit the sizes of the cen¬ 
tral bins, such that they are smaller or comparable to the 
seeing disc, but large enough to provide sufficient signal for 
a robust extraction of the kinematic features. 


4.2 MUSE kinematics 

We use the publicly available penalised P ix el Fit ting 
(pPXF)0 method of ICappellari fc EmsellemI ll2004ll to 
determine the LOSVD of stars. This method uses 

^ http://purl.org/cappellari/idl/ 
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Figure 1. MUSE spectra of the central region of NGC 5813, integrated within a circular aperture of a radius of 10^^, showing the pPXF 
(top) and GANDALF (bottom) fits. The bottom spectrum is shortened and shows only the fitted region. The fit is limited to the region 
blue ward of 7000A, as the spectra are significantly influenced by the tellurics and imperfect sky subtraction in the red. For the pPXF fit, 
we masked regions around the prominent emission (and sky) lines, but they are fitted with GANDALF. Legend shows the components of 
the fit and the prominent lines of interest. The pPXF fit gives a S/rN of 200, while GANDALF fit provides a S/rN of 255. 


the Gauss-Hermite p a ramet ris ation o f the LOSVD 
llvan der Marel fc Frarud 1 19931 : iGerhardl Il993l ). which 
are used to estimate departures in shape of the LOSVD 
from a Gaussian, such as the skewness and kurtosis. pPXF 
builds an optimal spectral template from a library of 
stellar (or gala:xy) spectra, which when broadened with a 
given Gauss-Hermite function, will best fit the observed 
spectru m. We used the entire MILE^^ library of stellar 
spect ra (ISanchez-Blazauez et al.ll2006l : iFalcon-Barroso et al.l 
I 2 OIII ). The resolution of the stellar library (2.50A, or cr = 48 
km/s at 6560A) is smaller, but very similar to the resolution 
of the MUSE spectra, which show a changing resolution: 


® MILES library is available from http://miles.iac.es 


at 5000 A of 2.74 (69 km/s), at 6000 A of 2.59 (54 km/s) 
and at 7000A of 2.54 (46 km/s). As the resolution is 
significantly lower than the stellar velocity dispersion in 
NGC 5813, we did not convolve the spectral library to the 
same resolution of the MUSE data. 

Both the library and the MUSE spectra were logarithmi¬ 
cally rebinned and the pPXE fit was performed in two steps. 
The purpose of the first step was to derive the optimal tem¬ 
plate by fitting the full spectral library (-^ 1000 spectra) to a 
representative spectrum, which was created by collapsing the 
full NGC 5813 data cube. We set up pPXF using 6 moments 
of the LOSVD (U, cr, h^, hi, / 15 , and he), and we added a 4th 
order additive polynomial. We used the default value for the 
penalisation (0.4). Finally, we limited the wavelength range 
used for the fit to A < 7000A, and excluded the regions with 
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Figure 2. Maps of LOSVD moments of NGC5813. Top row: the mean velocity, the velocity dispersion, Gauss-Hermite coefficient ha. 
Bottom row: Gauss-Hermite coefficients / 14 , /15 and fig- Colour scale is given by the colour bar at the bottom, while the ranges for each 
moment are shown on the lower-right corner of each map. When plotting the and hg maps we subtracted the medians of 0.012 and 
-0.015, respectively. The black contours are isophotes from the white light MUSE image. The straight line on the mean velocity map 
shows the orientation of the KDC and the inner photometric major axis (PA=146°). The dashed line on the velocity dispersion map 
shows the outer (global) photometric major axis (PA=135°), which coincides with the orientation of the 2a peaks. The diamond symbol 
is the location of the bin used in Section 14.31 and Fig.|4] North is up and East to the left. 


possible emission lines. We also excluded the region around 
the strong 5577 A sky line which was not removed well for all 
spectra. A fit to a spectrum extracted from a circular aper¬ 
ture is shown on the top panel of Fig. [T] In order to verify 
the extracted kinematics, we fitted various spectral regions. 
For example, we fitted the full wavelength range and the 
SAURON wavelength range (4570 - 5200 A) for comparison. 
These fits resulted In the same general trends and kinematic 
features. 

In the next step we run pPXF on the individual spectra 
of the MUSE data cube using always the optimal template 
constructed from the five spectra with non-zero weights as¬ 
signed by the pPXF fit in the first step. For each MUSE 
spectrum we estimated S/rN, where the rN (residual noise) 
was estimated as the resistant standard deviation of the dif¬ 
ference between the spectrum and the pPXF best fit, limited 
only to the wavelength region used in the fit (also exclud¬ 
ing emission-lines). The central regions are characterised by 
S/rN of around 130, which then steadily drops to about 90 
at the distance of 28”. 


The uncertainties were estimated running Monte Carlo 
simulations for each bin, where the original spectrum was 
randomly perturbed 500 times. The level of perturbation is 
determined from the rN value of the original pPXF fit to the 
spectrum. The only change to the pPXF set up was in the 
removal of the penalisation during the Monte Carlo simula¬ 
tions. Average uncertainties on V and a are 3 km/s and 4 
km/s, respectively, while for the /13 to he Gauss-Hermite mo¬ 
ments they are ~ 0.013. Errors are constant across the field 
with a dispersion of about 0.5 km/s for V and a, and 0.0025 
for the Gauss-Hermite moments, with a small tendency for 
an increase in the outermost (edge of the field) bins. 

Figure presents the kinematic maps of NGC 5813: 
the mean velocity, velocity dispersion and the higher order 
Gauss-Hermite coefficients: he —he. The M USE data confirm 
what has been observed a lready before dEfstathiou et all 
Il982l : lEmsellem et al.ll2004h : the KDC In the central 10” of 
the velocity map, followed by a zero level mean motion, and 
the complex structure in the velocity dispersion maps con¬ 
sisting of a central peak, a nearly circular region with a signif- 
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Figure 3. Comparison of MUSE (black circles) and SAURON 
(red open squares) velocity and velocity dispersion values along 
the photometric major axis (at large scales). Note the similar¬ 
ity of the two kinematics measurements, but the higher quality 
of MUSE data, both in terms of smoother variations and higher 
spatial resolution. 

icantly lower a values around it, the subsequent rise and the 
final decrease of the velocity dispersion values. The map 
shows anti-correlation with the KDC velocity structures and, 
as expected, the /15 map shows a correlation. The median val¬ 
ues of hz and /is maps are not zero, but 0.012 and -0.015, 
respectively. This is possible due to a minor template mis¬ 
match, but within the errors, the maps are consistent with 
zero outside the KDC region. For presentation purposes we 
removed the median values when plotting. /14 has a relative 
drop in the very centre (associated with the central rise in 
(j), surrounded by an elevated region corresponding to the 
KDC area, which then falls off towards the edges. The he 
map is essentially fiat. 

4.3 Kinematic evidence against two thin 
counter-rotating discs 

In the case of NGC 5813, the most striking is the structure on 
the velocity dispersion map. In Fig. [3] we show the similari¬ 
ties between the SAURON and MUSE velocities and velocity 
dispersions along the major axis of NGC 5813. The existence 
of the complex structure on the velocity di spersion maps is 
alrea dy clear on the SAURON kinematics llEmsellem et al.l 
|2004| ). but due to the shape of the field-of-view of the final 
SAURON mosaic, and the noise in the velocity dispersion 
map, it was not possible to recognise the specific shape of 
the high a structure. The MUSE data show that the region 
with the elevated velocity dispersion is symmetric and flat¬ 
tened along the photometric major axis (and not, for exam¬ 
ple, distributed in a circular ring). Furthermore, the velocity 
dispersion map consists of two distinct peaks along the pho¬ 
tometric major axis of the galaxy (the velocity dispersion 
is lower along the minor axis). The peaks occur around 11” 
from the centre, coinciding with the edge of the KDC, as can 
be seen in Fig. |31 

This shape is in a striking resemblance to the shapes 
of the velocity dispersion maps found in 2a galaxies. 
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Figure 4. Gaussian LOSVDs extracted at position (-7.5,-8.0) and 
10.9” away from the nucleus in the south-east region of the high 
velocity dispersion are. All LOSVDs are normalised such that 
their integral is equal to one. The single component LOSVD is 
shown by black full line, while the two separate LOSVD, per¬ 
taining to the prograde and retrograde motion are shown with 
blue (dashed) and red (dot-dashed) lines, respectively. A com¬ 
bined LOSVD (adding prograde and retrograde components) is 
shown by a green (dotted) line. The systemic velocity from a sin¬ 
gle component fit was removed from all LOSVDs. 

made of two counter-rotating discs (for examples of veloc- 
ity and velocity disp ersion maps see figs C.5 and C.6 of 
iKrainovic et al.l201ll ). As discs counter-rotate, the combined 
LOSVD is broadened and deformed, where the relative speed 
of rotation might even induce a clearly separated two peak 
shap e of the LOSVD (e.g. as in NGC 4550, see iRix et al.l 
I1992I ). In general, the broadening does not have to be suffi¬ 
cient for such a separation in two peaks (the relative speeds 
are not sufficiently different for such a separation), but in¬ 
stead manifests itself in a general broadening of the LOSVD 
along the major axis (e.g. as in NGC 4473). The suppres¬ 
sion of the mean rotation is also a direct consequence of the 
counter-rotation. If the two components have similar masses 
and distributions (e.g. disc thickness and scale lengths) then 
there will be no net rotation. If, however, one component is 
dominant within a certain region, the result will be an in¬ 
terchange of rotatio n, which can result in a counter-rotating 
core (e.g. NGC0448 IKrainovic et al1l201ll ). or a disappear¬ 
ance of rotation i n the central region (NGC 4550) or at larg er 
radii (NCG 4473. lEmsellem et al.ll20oJ : iFoster et al.ll^Ola) 
NGC 5813 is, however, probably not made of two 
counter-rotating, relatively thin, discs like a typical 2a 
galaxy. A typical 2a galaxy has ellipticity larger than 0.5, 
althoug h there are cases with 0 .3-0.4 ellipticites in the outer 
regions (|Krainovic et al.ll201ll ). NGC5813 shows a signifi¬ 
cant change in ellipticity with radius, being 0.1 within the 
KDC region and somewhat flatter at larger scales (ellipticity 
of 0.27). Furthermore, as seen on Fig.O the KDC is oriented 
at approximately 145°, but the two peaks in the velocity dis¬ 
persion are oriented at approximately 135°. These two orien¬ 
tations are also reflected in the varying position angle of the 
major axis, which changes from ab out 145° in the central 5 ” 
to 135° beyond the central 10” (e.g. lPullo fc Grahamll20l3 ) . 
It is certainly telling that the orientations of the kinematic 
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features are linked to those of the photometric body, but 
they are not compatible with two thin discs. Even if discs 
are misaligned (one at 135° and the other at 145°), it would 
be difficult to account for such significant cancellation of the 
observed mean motion outside of the KDC. If there are two 
counter-rotating components, they are likely not fast rotat¬ 
ing discs. 

If the LOSVDs of NGC5813 are indeed made of two 
counter rotating components, it could be possible to extract 
each co mponent separately, as for N GC4550 or some other 
galaxies dCoccato et al.|[201ll . I20f3l . In the cases where the 
disentangling of the components is successful, there is typi¬ 
cally a substantial difference in stellar populations, and, even 
if co-rotating, a substantial difference in kinematics (e.g. 
Ijohnston et al.l[20131 : iFhbricius et al.l[20l3 l. In this case, as 
is shown below, there seems to be no age difference between 
the prograde and retrograde kinematics components, while 
there is a certain difference in the metallicity. This makes the 
disentangling of the components more difficult, and therefore 
we attempted to separate the components only at one rep¬ 
resentative spatial location. 

We used the pPXF ability for two component fitting 
and, to probe the stellar population differences, we made use 
of the single st ellar population spect ra from the MIUSCAT 
model library (IVazdekis et al.l l2012l l , designated as safc0, 
spanning a range of ages 0.06 - 14 Gyr, and metallicities 
[Z/H]=-0.4 - 0.22. We first ran pPXF in a single compo¬ 
nent mode, using the full MIUSCAT library (200 spectra), 
and determined which stellar templates were used in the fit. 
Our selected spectrum can be fitted with four MIUSCAT 
templates, which are all old (11-14 Gyr) and have a typi¬ 
cal metallicity of either -0.4 or 0.2. From these we selected 
two with significant weights, and the largest difference in 
the metallicity (ages did not differ), assigned then to each 
component and ran pPXF in the two component mode. In 
all cases we used Gaussian LOSVDs. We selected a bin in 
the south-east region of the high velocity dispersion values 
(marked by a diamond in Fig. [2]). 

Resulting LOSVDs are shown in Fig. |4] There are: the 
LOSVD of the single component fit (similar to the kine¬ 
matics presented in Fig. [21 but limited to only the first two 
moments), retrograde and prograde LOSVDs from the two 
component pPXF fitfl, and a combined LOSVD of the retro¬ 
grade and prograde components. As one would expect, the 
single component LOSVD is rather similar to a combined 
LOSVD of the retrograde and prograde components, and 
this explains the characteristics of the observed V and a 
maps. The two components are different in their velocity 
dispersions and mean velocities. The difference in the lat¬ 
ter, for this spatial bin, is only about 100 km/s, and much 
less than the velocity dispersions (326 km/s and 196 km/s 
for prograde and retrograde components, respectively). The 
large velocity dispersions of both components and their rel¬ 
ative low separation in the velocity indicate that NGC5813 
is not made of two discs. We postpone the discussion on the 
implication of this result to Section [T] but stress that the 
process of spectral disentangling is highly degenerate, espe- 


® MIUSCAT library is available from http://miles.iac.es 
^ We define the retrograde component as the one that has a sys¬ 
temic velocity lower than the single component fit. 


daily If there is no significant difference in stellar population, 
or the separation in velocities is small. Therefore, the results 
presented in Fig. U should be interpreted with caution. 


5 AN APPROXIMATE DYNAMICAL MODEL 

OF NGC5813 

If NGC 5813 kinematics are a result of two counter-rotating 
structures, orbit-based dynamical models can be used to de¬ 
duce the relative mass fractions of the components and the 
distribution of the orbits that constitute them. We use the 
ISchwarzschildl (| 19791 . 0*98211 orbital superposition method for 
constructing asymmetr ic ga laxies, with the code de fined in 
ICappellari et al] (l2006l l and ICappellari et al.l ll2007h . Those 
papers explain the workings of the code in detail, which we 
do not repeat here, but give a brief overview of the necessary 
steps and choices for starting parameters. 

5.1 Mass models and the Schwarzschild three 
integral method 

The essence of the Schwarzschild method is the construction 
of a library of orbits, possible in the gravitational poten¬ 
tial defined by a mass model. The orbits from the library 
are then combined building a model galaxy, which is pro¬ 
jected on the sky, calculating the observables comprising 
the mass distribution and the histograms of the LOSVDs at 
each position on the sky. These are compared to the original 
mass model and MUSE kinematic observations. The model 
is constrained using the non-negative least-square method 
llLawson k, Hansmilll974l L Evidently, there are three crucial 
ingredients needed to robustly constrain the dynamical mod¬ 
els: a mass model, a representative orbit library, and high- 
quality kinematics. 

We use the parametrisat i on of fight of NGC 5813 from 
Table B2 of ICappellari et ^ (l2006l L which is based on the 
Multi- Gaussian Expansion method (MGE, lEmsellem et al.l 
11994311 and constructed from an HST/WFPC2/E814W im¬ 
age and a ground-based image from the MDM observa¬ 
tory in the same filter. The step from the MGE to a mass 
model is accomplished by deprojectin g the observed surfac e 
brightness using the formulas given inlhfonnet et al.l (Il992l l , 
lEn^ellem et al.T(^1994b^ .lCappellaril ll2002ll. as well as the ro¬ 
bust approximations from Cappellari et al.| (l2002l l , assuming 
axisymmetry and a certain (constant) mass-to-light ratio, 
M/L. The stellar orbits are then calculated within the po¬ 
tential defined by the stellar density and a central massive 
black hole. For the deprojection we assume the inclination of 
89°, but we also run a model with inclination of 60°, which 
gave essentially the same results. 

The orbit library needs to be sufficiently large to pro¬ 
vide a complete set of brisis function that can describe the 
distribution function of a galaxy. The fundamental problem 
is that the orbits are essent ially ^-functions while t he distri¬ 
bution function is smooth. ICappellari et al.l (l2006l l resolved 
this issue by adding up a large number of ’dithered’ orbits 
started from adjacent initial conditions. Each orbit is defined 
with three integrals of motion: energy E, the ^-component 
of the angular momentum Lz and the third integral I 3 . The 
orbits are here sampled on a three dimensional grid with 
41 energy points (corresponding to the representative radius 
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of the orbit), and at every energy we used 10 angular and 
10 radial sectors, arranged in a polar grid, with the number 
of ’dithers’ equal to six for each dimension of the grid (for 
detailed visualisation of the or bital starting scheme an d the 
dithering pattern see fig. 6 of ICappellari et al.1 l2006ll . Ac¬ 
counting for prograde and retrograde orbits, the complete 
orbit library consist of more than 1.7 million orbits. 

The potential within which the orbits are integrated 
is defined by the choice of symmetry, mass model and ad¬ 
ditional (dark) sources of gravity (central black hole and 
the dark halo). We run the axisymmetric version of the 
Schwarzschild code limi ting the possib le orbits to various 
sorts of short axis tubes (Ide Zeeuwlll985l l, but other resonant 
and chaotic orbits can occur. The observations presented in 
Sectionsl^and fOl indicate that the galaxy is at least weakly 
triaxial. Our approximation is, however, sufficient for the 
purposes of the present study, as we are not interested in the 
detailed shape of the galaxy. We want to test if the MUSE 
kinematics (the KDC and the 2 (t peaks structure on the 
velocity and velocity dispersion maps, respectively) can be 
reproduced by two counter-rotating orbital families, which, 
given the orientation of the KDC, have to rotate around the 
minor (short) axis. As much as they could be present, the 
long-axis tubes and box orbits are not expected to contribute 
much to the observed kinematics, and we are not interested 
in determining the mass of the black hole (M;,^), dark mat¬ 
ter contribution to the overall M/L, or the three dimensional 
shape of the galaxy. 

Finally, crucial constraints for a Schwarzschild model 
are provided by the MUSE kinematics (Section 3)). As the 
model is by construction axisymmetric, we symmetries the 
MUSE data using point-(anti)symmetry around the global 
photometric axis (PA=135°), averaging the values at posi¬ 
tions: [(x,y), (x,-y), (-x,y), (-x,-y)], but keeping the original 
errors for each point. 


5.2 Orbital space of NGC 5813 


We construct a Schwa rzschild model of NGC 5813 assum¬ 
ing M/L=4.81 (as in ICappellari et al.] [2Q(M) and a M^^ 
from the recent Mbh - u relation of McConnell &: m3 
20lj). Using the effect ive velocity dispersion of 210 km/s 


I Cappellari et alll2013bll . the expected Mi,^ = 2 x 10® Mq. 


At the distance of 31.3 Mpc this black hole will have a radius 
of the sphere of influence of less than Cf.'2, which is below our 
resolution elements (the seeing was about O'.'? FWHM). Al¬ 
though Schwarzschild models constrained with high quality 
IFU data can probe s cales about 3 times sm aller than the 
seeing resolution (e.g. iKrainovic et al.ll2009l l. we do not fit 
for the Mbh- Therefore, we are not concerned about the ex¬ 
act shape of the PSF, which we parameterise with a single 
Gaussian (o-psf = 0.3"). We are primarily interested in the 
region between about 1-30", the region dominated by the 
KDC and the outer, non-rotating parts of the galaxy which 
show the 2a structure. We also ignore the contribution of 
the dark matter halo. The half-light radius of NGC 5813 is 
about 57", and the MUSE field-of-view, as large as it is, 
still samples just above a half of the effective radius. The 
contribution of the dark matter within t his region is ex¬ 
pected to be moderate (~15 per cent), but ICappellari et al.l 
ll2013al l find about 47 per cent of dark matter within th e 
half radius. Their Jeans anisotropic model (ICappellarill2008| l . 


however, has a quality flag 0, indic ating a problematic data 
or model llCappellari et aDl2013bl l. The comparison of the 
model, which reproduces well all major kinematics features, 
with the data is shown in Fig. I All of Appendix [Al 

In Fig. [5] we show the orbital space of the Schwarzschild 
model. There are 19 panels, each showing a polar grid of 
orbital starting position for a given energy (represented with 
the radius of a circular orbit at that energy). The radii were 
chosen to sample the KDC and the surrounding region. Each 
panel is divided in the right and left side, corresponding to 
prograde and retrograde orbits, respectively. The definition 
of what is prograde and retrograde is arbitrary, and we assign 
the component responsible for the KDC to be the prograde 
component. High angular momentum orbits are found in the 
lower right and left corners of the panels. Orbits with an 
increasing /a values are found higher above the x-axis. The 
colour contours show the relative mass weights assigned to 
orbits in the construction of the model. 

Two important results can be seen in this figure. Firstly, 
in every panel there is a significant fraction of mass assigned 
into both prograde and retrograde orbits, but this fraction 
is changing with radius. Within the region of the KDC, and 
especially the region where the observed rotation curve is 
rising (<5"), the orbital space is dominated by the prograde 
orbits (~ 70 per cent). This ratio diminishes towards the 
visible end of the KDC, where about half of the orbits have 
one and the other half the opposite spin. In the region of 
the 2cr peaks (~10 - 25"), and at the largest radii probed 
by the model (~30"), the retrograde component is slightly 
more dominant. 

Secondly, the orbital space of NGC 5813 lacks high an¬ 
gular momentum orbits, located at the right and left corners 
of the semicircles. Instead, at each energy, orbits are defined 
with both moderate and I 3 integrals. There is, however, a 
difference between the two components, where the retrograde 
orbits are characterised with a typically smaller angular mo¬ 
mentum compared to the prograde orbits. In all panels the 
retrograde orbits are confined to the region with R/Rc < 0.5, 
while this is not the case for the prograde orbits. This implies 
that the prograde component is overall rotating faster. This 
is also likely what drives the slight Imbalance in the mass 
fractions beyond the KDC, where the net rotation needs to 
be zero. 

We run several models with different starting param¬ 
eters: the inclination (also at 60°), the orientation of the 
symmetry axis (also at PA=145°, the orientation of the in¬ 
ner photometric major axis and the KDC), the black hole (an 
order of magnitude smaller and higher mass) and the M/L 
ratio (with A M/L=l). We also constrained the models with 
the kinematics extracted at the target S/N of 70 (same as 
used for the ionised gas extractions). All these changes gave 
rather similar orbital configurations always characterised by 
two components, where the ratio of the prograde to retro¬ 
grade orbits changes from the initial 70:30 in the centre to 
50:50 towards the end of the KDC region, while in the outer 
regions the orbital families had a 45:55 ratio. 

This implies that the best fitting orbital space is made 
of two, relatively thick, short-axis tube orbit families, which 
permeate the system. The structures are counter-rotating, 
but they are not thin discs made of stars at high angular mo¬ 
mentum orbits with low third integral. Instead they are co- 
spatial vertically thick components, approximately identical 
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Figure 5. Integral space of an axisymmetric orbit-superposition model for NGC5813. Each panel plots the meridional plane {R,z) for 
a given energy with the starting positions of orbits shown with white dots. Energies are presented with the radius of the corresponding 
circular orbit in arcsec (green numbers in the upper left corner). Prograde orbits are shown on the right side (with a ‘-I-’ sign), while the 
retrograde on the left side (with a ’ sign) of each panel. The radius of the circular orbit is the size of the horizontal axis (measured 
from the centre). Orbits with the highest angular momentum are found in the right and left corners. Over-plotted is the fraction of mass 
assigned to orbits at constant energy. The percentage of the total mass assigned to prograde and retrograde orbits are printed above the 
corresponding side of the panel. We show only a subset of radii, covering the KDC and the outer body fully within the MUSE observations. 
The colour bar represents the relative colour coding of the orbital mass weights, where 1 represents the largest mass weight assigned to 
an orbit at the given energy. 


in mass in the outer parts, but with the prograde component 
dominating in the centre (within the KDC). In Section [3 we 
will discuss further the implication of this orbital structure. 


6 EMISSION-LINE KINEMATICS 

NGC 5813 is known to have an extended a nd morpho¬ 
logically complex distrib ution of ionised gas llSarzi et al.l 
I2OO6I : [Randall et al.ll20lil ). and it has a complex filamentary 
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Figure 6 . MUSE ’narrow’ band image of Ho emission in 
NGC 5813. The image as obtained by summing the flux within 
a narrow range around Ho + [Nii] complex (between 6615 — 6640 
A), estimating the continuum level around that region and di¬ 
viding it out. North is up and East to the left. Contours are the 
isophotes from the MUSE white light image. 

dust distribution with an apparent nuclear dust disk (e.e . 
ICaroIlo et al.|[l997bl : iTran et al.ll200lh . IWerner et al.l (l2014l ( 
detect also cold (100 K) [Cll] gas, which is co-spatial with 
the Ho emission-line as well as the lowest entropy X-ray hot 
gas. The complex and multi-phase gas content of NGC 5813 
is most likely re lated to the cooling a ssociated with the AGN 
activity cycles llRandall et al.ll201ll i. We investigate here if 
there is a connection between the gas content and the stel¬ 
lar KDC. A more thorough investigation of the emission-line 
gas properties we postpone to a future study. 

6.1 Distribution and kinematics of emission-line 
gas 

In Fig.l^we show the Ha + [Nil] image summing the MUSE 
data cube between 6615 and 6640 A. The strong Ha + [Nil] 
features consist of two major filaments, one in the North- 
East direction, and the other in the South-South-West di¬ 
rection. Some other, less pronounced, filaments and knots 
are also visible, notably in the East part of the map. Other 
emission-lines typically follow the features seen here in Ha 
+ [Nil] emission, although some are less prominent. Having 
the spectral information, we can proceed further and mea¬ 
sure the kinematics o f ionised ga s. _ 

W e use Gandalf llSarzi et al.ll200d : lEalc6n-Barroso et al.l 
I2OO6II to measure fluxes and kinematics of the fol¬ 
lowing emission-lines: H/3, [Oiii[AA4959,5007, [Ni[AA5198, 
[Ni[AA5200, [Oi[AA6300,6364, Ha, [Nii[AA6547,6583, and 
[Sll[AA6716,6731. Eig. |B]shows that emission is abundant in 
certain regions and therefore we decided to change the bin¬ 
ning from the high S/N used for the st ellar kinematics. We 
again used the Voronoi-binning code of ICappellari &: CopinI 
(I2OO3I I setting the target S/N to 70 (again at 5500 A). This 
S/N is sufficient to robustly estimate the absorption features 
before attempting to fit the emission-lines, but it results in 
smaller bins: now we have more than 4000 spectra, which 


within the central 10" are rarely binned and have the orig¬ 
inal 0'.'2 sampling. The central S/rN (of the pPXF fit) to 
these spectra approaches 140. 

The extraction of emission-line kinematics consisted 
first of a fit to the stellar continuum using pPXF and 
the MILES library, followed by a fit to the observed spec¬ 
trum in which the optimal template (obtained from the first 
pPXF ht) is combined with several Gaussians mimicking the 
emission-lines. The call to pPXF differed somewhat from the 
one in Section S31 as we parameterised the LOSVD with 
only four moments and used the 10th degree multiplicative 
polynomial. The fitted region was the same as for the kine¬ 
matics. An example of a GANDALF fit to a spectrum within 
a central circular aperture is shown on the lower panel of 

Fig.Il 

We set up Gandalf to fit freely only Ha and the stronger 
line of the [Nil] doublet, while the kinematics of other for¬ 
bidden lines were tied to the [Nil] line and of H/3 to Ha. 
This means that kinematic information could potentially be 
different only for Ha and [Nil] doublet, but fluxes of all lines 
were independently estimated. 

We show the distribution of all detected lines in the 
Appendix iB] (Fig. IBIII . while in Fig. [7] we show the fluxes, 
velocities and velocity dispersion for Ha and [Nil]. Only 
those bins for which the lines had an Amplitude-to-residual- 
Noise (A/rN) ra ticfl greater than 3 are plotted (following 
ISarzi et al.ll200d L With this criterion all of the fitted lines, 
except [Nl], were detected, but [Ol] was detected only in a 
very limited region, still following the general distribution 
of the [Nil] line. The distribution of Ha or [Nil] on Fig. [3 
closely resembles what is seen on the reconstructed image 
in Fig. |S] Similar structure was also seen in the SAURON 
maps of [Olll] an d H/3 emission-lin es, but at ~ 5 times lower 
spatial sampling dSarzi et al.ll200^L 

Ha and [Nil] lines show similar and complex kinematic 
structures. Compared to the systemic velocity of the galaxy, 
the ionised gas in the southern arm is receding, while the 
northern arm shows both receding (east side of the arm) and 
approaching (west side of the arm) material. The measured 
velocity dispersion is relatively uniform over both arms and 
other filaments, typically between 60 km/s and 120 km/s 
(note that the instrumental resolution at the wavelength of 
Ha is about 50 km/s and was not removed from the maps). 
The highest velocity dispersion values, in the excess of 220 
km/s are found slightly off-centred towards the west. 

Emission-line gas is distributed almost in a polar con¬ 
figuration, extending close to the minor axis of the galaxy, 
and co-spatial with direction of the radio emission and 
the location of buoy ant cavities in the X-ray emitting gas 
([Randall et al.|[201lh . The distribution of the emission-line 
gas in the southern region is marked by two regions with no 
detections. These cavities in the emission-line gas are cen¬ 
tred on approximately (5",-5") and (15", -15"), and they are 
co-spatial with the two innermost cavities in the X-ray gas. 
Given that the cold gas is most likely p roduced by thermally 
unstable cooling from the hot phase, IWerner et al.l (l2014h 
suggest that the filamentary emission-line gas is linked to 
the buoyant cavities and the jet activity, where the relativis- 


® The amplitude refers to the amplitude of the fitted Gaussian. 
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Figure 7. Emission-lines Ho (top) and |Nii] (bottom) in NGC 5813. From left to right: A/rN map where only emission-lines with A/N> 3 
are shown (region with no reliable emission is coloured grey), logarithm of emission-line flux (in units of 10“^® erg/s/cm^), the mean 
velocity and the velocity dispersion of a single Gaussian. Contours are the isophotes from the MUSE white light image. 


tic plasma interact with the filaments and drags them out 
from the centre. 


6.2 Central emission-lines kinematics 


Both Ho; and [Nil] emission-lines show peaks in the velocity 
dispersion in the very central region, not coinciding with the 
centre of stellar distribution (which is about 1” to the West). 
The peaks in emission-line velocity dispersion also coincide 
with a region which shows a different velocity structure of the 
ionised gas compared to the outer parts. Zooming in onto the 
central 10"xl0”, the different kinematics in the central 3” 
becomes apparent. We show this on Fig. (HI Coinciding wit h 
the nuclear dust disc seen on HST images I Tran et al.ll200l] ). 
somewhat more extended (~ 3”), but with nearly the same 
orientation, there is a kinematic structure with approaching 
side in the North-West and the receding side in the South- 
East. The approaching North-West side seems limited to the 
central 2”, surrounded by receding material. This interface 
between the approaching and receding material is also char¬ 
acterised by the peak in the velocity dispersion. The South- 
East side is more extended but also not uniform, showing 
two regions of increcised velocities. This complex is present 
in both Ha and [Nil] emission-lines and displays similar dis¬ 
tribution and kinematics. Given that it is spatially associated 


with the dust disc, it is possible that the velocity structure 
has the same inclination as the dust disc (nearly edge-on). 
The disc nature of the feature might, however, be difficult to 
reconcile with the off-centred high velocity dispersion peak 
and the non-regular spatial distribution of velocities men¬ 
tioned above. It is likely that the observed motions are not 
solely of gravitational origin. 

The higher S/N of MUSE data and the small size of the 
structure are possibly the reasons why this kinem atic feature 
was n ot seen previously in the SAURON data llSarzi et al.l 
l200d) . Alternatively, it could also be a feature that is as¬ 
sociated more with Ha and [Nil] than with the lines ob¬ 
served within the short SAURON wavelength range (e.g. 
[Olll] which was used to derive the kinematics). In the MUSE 
data, H/3 and [Olll] are essentially not detected in the very 
centre, although there is a hint of the kinematic structure 
on the [Olll] maps (see maps on Fi g. IBlIl . This might also 
explain whv lMcDermid et al.l (|200(il ) . who presented higher 
spatial resolution observations of SAURON sample galaxies 
with OASIS, also did not see this feature. 

Of particular interest is the fact that the central ionised 
gas rotation is oriented similar to the KDC and the nuclear 
dust disc (shown with a line on Fig. [8]). We stress that the 
region of the high velocity dispersion in both Ha and [Nil] 
emission does not coincide with the centre of rotation, but 
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Figure 8 . Kinematics of Ho (top) and [Nil] (bottom) emission¬ 
lines within the central 10 X 10”. The mean velocity is on the left 
and the velocity dispersions on the right. The solid line on the 
velocity maps indicates the coinciding orientations of the stellar 
KDC and the dust disc visible on the HST images. The dashed 
vertical and horizontal lines indicate the position of the galaxy 
centre. Contours are the isophotes from the MUSE white light 
image. 


If this scenario is correct, we are not seeing an inflow 
of the gas towards the centre, which settles in a disc-like 
structure and eventually creates stars and contributes to the 
stellar content and the kinematics of the KDC. On the con¬ 
trary, we are seeing a gas reservoir which resides in the equa¬ 
torial plane of the galajcy, and which could originate from the 
stellar mass loss, a relatively recent accretion or even be a 
remnant of a gas disc, which was, at an early epoch, re¬ 
sponsible for the creation of the prograde rotating structure 
and the KDC. Currently this gas could be serving as the 
AGN fuel reservoir. The jet and the buoyant X-ray cavities 
are lifting and dragging the material from this gas reservoir, 
which is being cooled and forms a multiphase gas structure, 
as observed in the cold [Cll] or warmer Ha and [Nil] emis¬ 
sion. The velocity structure of the northern Ha arm actually 
connects to the central gas disc, possibly supporting the ev¬ 
idence that gas is being lifted from the surface of the disc, 
but retaining its angular momentum. 

Therefore, the observations suggest that the present gas 
in NGC 5813 is not of high relevance for the KDC. Further in¬ 
vestigations of the properties of the emission-line gas, which 
could be used to asses this scenario, are beyond the scope of 
this paper, but we note that high spatial resolution observa¬ 
tion of the cold gas phase centred on the nucleus would be 
beneficial. 


7 DISCUSSION 

Previous sections presented MUSE observations of 
NGC 5813 and their analysis. Here we summarise the 
evidence for the counter-rotating nature of the KDC and 
discuss the possible formation scenarios of NGC 5813. 




surrounds the region of the approaching side of the central 
kinematic feature. As it arises where there is a reversal in 
the velocities, this could indicate shocks and the turbulent 
nature of the interstellar medium, linked with the relativistic 
plasma of the jets and the buoyant cavities. 

6.3 The multiphase gas, the activity and the KDC 
of NGC 5813 

The central, possibly disc-like, kinematic structure revealed 
by MUSE is different from the rest of the ionised gas dis¬ 
tributions. Its orientation is perpendicular to the radio jet 
and the direction of rising X-ray cavities, but it is oriented 
similar like the KDC. Is this gas feature related to the KDC? 

Stars within the KDC region are old and there is no 
evide nce for any significant c urrent star formation in this re¬ 
gion (iKuntschner et al1l2010l ). The cold gas that is observed 
in NGC 5813 is likely a result of cooling instabilities, which 
also seem to be typ ically present in cold-gas-rich systems 
(IWerner et al.ll2014l ). The spatial link between the orienta¬ 
tion of the jet and the cavities filled with relativistic particles 
and the emission-lines gas detected in MUSE observations, 
suggests that these gas filaments are being dra gged outwards 
by th e jet and the bubbles. As described in I Werner et al.l 
(|2014| ). the hot particles can penetrate the filaments and 
heat the cold gas, both contributing to its distraction and 
creation, depending on their energy. 


7.1 The counter-rotating nature of the KDC in 
NCG5813 

The stellar kinematics of NGC 5813 observed with MUSE 
suggests that this galaxy is made of two counter-rotating 
components. This is primarily based on the velocity disper¬ 
sion map, which features an elongated structure with two 
peaks along the major axis of the galaxy, located some 11” 
from the nucleus. The appearance of the u-peaks coincides 
with the visual end of the rotation and the KDC. In their 
spatial distributions, the velocity and the velocity disper¬ 
sion maps are similar to maps of 2 a galaxies, for which it 
has been shown that they are made of two counter rotating 
structures. 

The stellar populat ions of NGC 5813 are old everywhere 
(iKuntschner et ciI]l201Cll ) . The pPXF fit based on stellar pop¬ 
ulation models selects typically a few old stellar templates 
for reproducing the observed spectra, which, however, differ 
in metallicity. Fitting two templates of old ages, but different 
metallicities to a spectrum in the region of the a peaks, re¬ 
turns two velocity components, with mean velocities on the 
opposite sides of the systemic velocity for that spectrum. 
While this method is strongly degenerate, it suggests that 
there are two physical components with different sense of 
rotation. 

An orbit-superposition model, which fits the stellar 
kinematics well, provides further evidence for the two compo¬ 
nent nature of NGC 5813. The two model components have 
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opposite spins, where the prograde component typically has 
a higher angular momentum on average, Finally, the mass 
fractions assigned to the orbits vary with radius, but both 
orbital families are present everywhere. Within the KDC re¬ 
gion, the prograde component is more dominant, explaining 
the KDC appearance, but in the outer regions the compo¬ 
nents contribute with essentially the same fraction of mass 
(small differences could be attributed to the difference in the 
angular momenta). 

The KDC region is also characterised by a drop in the 
velocity dispersion and a strong /13 moment, which is anti¬ 
correlated with the mean velocity. This is typically seen 
in systems where the LOSVD can be constructed of two 
components, a fast and a slow rotating ll Bender et al.|[l9^ : 
IScorza fc Bender! 1 1995lf . both in cases of co- and counter¬ 
rotation. This is consistent with the orbital space of the 
Schwarzschild model, which shows that the prograde com¬ 
ponent within the KDC has a higher angular momentum 
than the retrograde component. 

While the evidence for two counter-rotating structures 
in NGC 5813 are strong, it is also quite clear that these struc¬ 
tures are not counter-rotating thin discs, like in a typical 2 a 
galaxy. NCG5813 is more massive than a typical 2a galaxy. 
The shape of NGC 5813 is also not as flat as in known 2a 
galaxies. Moreover, the best htting orbital space is not popu¬ 
lated by high angular momentum orbits. Therefore, in spite 
of the fact that the outer light profile is best fit with an 
exponential (or almost exponential), the observed kinemat¬ 
ics are not well described by a large-scale disc, or discs. No 
net rotation in the outer regions would be compatible with 
a disc only if it is seen exactly at an inclination of zero de¬ 
grees (face-on). To show rotation, the KDC would then need 
to be made of a distinct component seen at a higher incli¬ 
nation. However, assuming the flattening within the KDC 
region (e ~ 0.1) is due to an inclined disc, the disc would be 
viewed at an inclination of about 25°, while the outer disc 
would be required to be at about 45° (e ~ 0.3). The stability 
of such a configuration even in a triaxial system is question¬ 
able. Furthermore, the two LOSVD components both have 
relatively high velocity dispersions, also arguing against the 
presence of thin discs. 

A simple solution is that the KDC is an orbital compos¬ 
ite, a result of a mixture of two components with different 
distribution of masses and similar, but not exactly the same, 
distribution of counter-rotating orbital families. A similar 
finding was already reporte d for another kinematically spec¬ 
tacular galaxy: NGC 4365 (Ivan den Bosch et al.|[200^ . This 
system is more complex than NGC 5813, as the axes of ro¬ 
tation of the KDC and the outer body are perpendicular to 
each other, where the outer body is in a prolate-like rota¬ 
tion. Still, the triaxial models of van den Bosch et al. show 
that stars that make the KDC in NGC 4365 are part of pro¬ 
grade short axis tubes which permeate the full system and 
the KDC is visible only due to an unbalance in the rela¬ 
tive weights in the central 6” between the various orbital 
families. 

The origin of the unbalance in the known orbital dis¬ 
tribution of KDC is not clear. A variety in possible for¬ 
mation processes for KDCs is highlighted by the range in 
the ste l lar populations ife ender &: Surmalll992l : ICarollo et al.l 
[1997^^ KuMschnerrtayj2010|) and kinematic properties 
dKrainovic et al. 20081 . 1201 ll ) of galaxies with KDCs. What 


then are the possible formation scenarios for NGC 5813 and 
its KDC? 


7.2 Assembly of NGC 5813 

KDCs were heralded as show cases of merger driven for¬ 
mation scenarios. Galaxies harbouring KDCs do not differ 
from other, similar mass, slow rotator early-type galaxies. 
Still, there are significant characteristics of NGC 5813 which 
point out to some possible formation pathways. 

The spatially uniform and old age of the stellar popu¬ 
lation in NGC 5813 argues for an early assembly, or alterna¬ 
tively, a later assembly of the system comprising star s cre¬ 
ated at high redshift {z ~ 3.5, Ku^schM^^^ 2^^. The 
negative metallicitv gradient j Efstath iou fc GorgasI 1 19851 : 
iGorgas et al.|[T99(]l : iKuntschner et al.ll201(j). constant abun¬ 
danc e of the a elements relative to iron (IKuntschner et ^ 
I 2 OIOI) . together with the finding in this paper that the 
two components permeate the (observed) body of the 
galaxy, argue against scenarios in whic h the KDG is a 
remnant of an accreted c ompact galaxy dKormen d^ 11984 
iFranx &: Illingworthlll98§) . 

Furthermore, kinematic observations presented here ar¬ 
gue against the origin of the KDC as a compact high-z el¬ 
liptical which grows in size (and mass) via minor mergers, a 
scenario invoked to explain the rapid size evolution of qu i- 
escent galaxies (e.g. Naab_et_aL 1200^; Ho2ldns_e£^jJb010|), 
but see als o (e.g . ICarollo et al.l 20131 : IPoggianti et al.ll20l4 
I Belli et aP I2OI4I) for an alternative scenario. A growth of 
NGC 5813 via intermediate to minor mergers (e.g. mass ra¬ 
tios of 1:5 to 1:10) is, however, not supported by the outer 
exponential profile of this galaxy. iHilz et aDd2013l ) show that 
such merger should increase the Sersic index to high values 
(n> 5), comparable to those found in giant ellipticals. 

Processes that could be considered in the case of 
NGC 5813 include: 


(i) Accretion of gas from a gas-rich companion which 
form s a disc-like structure in a principle plane of the system 
fe.g. iFranx &: Illingworthlll988ll . This pathway is invoked for 
the creation of 2a galaxies. An example is presented by 
ICoccato et al.l d201ll) . who show evidence for two counter¬ 
rotating discs in NGC 5719, a g alaxy which is stripp ing gas 
from its companion NGC 5713 dVergani et ahllioOTll . How¬ 
ever, there is no strong kinematic evidence for thin discs 
in NGC 5813, and NGC 5813 is a more massive system. 
Furthermore, the Schwarzschild model suggests a lack of 
high angular momentum orbits, and an existence of orbits 
with high vertical extent (above the equatorial plane). A vi¬ 
able pathway, however, is an early formation of the counter¬ 
rotating thin disc created from accreted gas, which due to 
subsequent interaction with satellites in the (sub)group envi¬ 
ronment was destroyed (dynamically heated and thickened). 

(ii) A major merger (mass ratios of 1:1) between two gas- 
rich spirals, or between a spiral and an early-t ype galaxy (e.g. 
Bender fc Surm3ll992l ) . Recent simulatio ns (iHoffman et al.l 


201(]l : Bois et al. j201ll : lMoodv et al.ll2014l) show that binary 


mergers of galaxies, where at least one progenitor has a mod¬ 
erate gas content (e.g. 15-50 per cent), can produce KDCs 
and 2cr galaxies. Subsequent m ergers of such gal axies can re¬ 
sult in disruption of the KDC ((Bqisj^,^ 201^, oreven an 
increased probability to host them (iMoodv et al.|[20l3 l . This 
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should be taken into account as NGC 5813, a central gala:xy 
in a (sub)group, likely experienced a few similar mass merg- 
ers. The analys i s of t he orbital structure of the remnants by 
iHoffman et al.l (l2010l ) reveals that KDCs are indeed a con¬ 
sequence of the balance of relative contributions of o rbital 
types spread throughout the system. iBois et ^ ll201lll show 
that if the spin of one of the progenitors is retrograde with 
respect to the orbital spin of the merger, the resulting KDC 
is also an orbital composite, with 2 a velocity dispersion fea¬ 
tures. Even initially prograde orientation of spins can result 
in a KDC, due to strong reactive forces which can cause the 
change in the orbit al spin directly before the final coalescence 
(iTsatsi et al.llioisl '). 

(iii) A formation as a natural consequence of the the hi¬ 
erarchical clustering scenario, where the galaxy sits on an 
inter section of filaments bringing cold gas l|Keres et al.ll2005l . 
I2OO9I I. A scenario was put forward as a case for the creation 
of counter-rotating disc galaxies (essentially the 2cr galaxies), 
whe re two streams are fe eding the host dark matter halo with 
gas dAlgorrv et al.|[2^014fl . If there is a sufficient difference in 
time between the arrival of the material (i.e. first through 
one filament and then through the other), two discs with op¬ 
posite spins can be created. A simultaneous feeding would 
likely result in a single component due to the collisional na¬ 
ture of gas, if the gas has not already fragmented into molec¬ 
ular cloud complexes. This scenario is somewhat problematic 
for NGC 5813 as there is no evidence for an age difference 
of stars. However, if the two feeding events happened early, 
were relatively short in duration and occurred within a few 
Gyr, it is essentially impossible to distinguish between the 
stellar populations. Moreover, given that NGC 5813 is one of 
the two most massive galaxies in the NGC 5846 group (and a 
centre of a sub-group) it is natural to imagine that the pro¬ 
genitor of the present day galaxy was formed on a location of 
merging cosmic filaments. The penetrating cold gas streams 
are rich in material and drive strong instabilities creating gi¬ 
ant s tar forming clumps within the galaxy disc (|Dekel et al.l 
I2OO9I I. Two streams, feeding the galaxy with material of op¬ 
posite angular momenta could produce two thick and turbu¬ 
lent discs, which through a subsequent evolution could start 
resembling the two orbital families seen in NGC 5813. 

All these scenarios require gas to be present during the 
formation event. The gas-rich formation scenarios have dif¬ 
ficulty accounting for a partially depleted core in NGC 5813 
nuclear surface brightness. The currently favourite scenario 
for creation of cores requires a bi nary black hole inter¬ 
action in a dissipation-less m erger dSegelman et al.l Il98fll : 
iMilosavlievic fc Merritt! 1200il l. The core might have been 
created later, during the subsequent evolution of NGC 5813 
as a central galaxy in a group. The constraint to this sce¬ 
nario is that the merger (or mergers) can destroy the origi¬ 
nal disc-like structure, but needs to keep the counter-rotation 
present. Furthermore, the existence of fast rot ators with par¬ 
tially depleted cores dKrainovic et al.ll2013bl') suggest there 
could be more unexplored pathways in creating kinemati¬ 
cally disc-like structures and removing stars from the nuclei. 

NGC 5813 is a massive galaxy, containing stars con¬ 
sistent with being created early, within a short period of 
time. As a central galaxy in a (sub)-group it s hould have 
expe rienced both phases of the mass assembly dOser et al.l 
l2010l l. In cosmological studies, galaxies with KDCs are still 


rare dNaab et al.ll2007l . l2014ll . found in remnants of systems 
that went though a dissipation-less similar mass merger. 
When analysed, these systems have various types of or¬ 
bits (short axis tubes, long axis tubes and box orbits), 
and the KDC is seen as an orbital composite dJesseit et al.l 
l2007l : iRotteers et al.|[20l3l . Slow rotator merger remnants, 
however, do not fully resemble the observed galaxies: they 
are typically too flat and, in particular, the velocity maps 
of KDCs are not yet comparable to the observed systems 
dNaab et al.l20f3l . More work is required here, and a promis¬ 
ing pathway in creating counter-rotating components (and 
KDCs in appearance) could be through cold streams which 
can create thick discs. The hope is that the halo mass and 
the properties of the streams will influence the outcomes, 
potentially creating a range of systems from less massive 2cr 
galaxies (with thin counter-rotating discs) to more massive, 
NGC 5813-type galaxies with KDCs (with thick counter¬ 
rotating discs). 

The case of NGC 5813 (together with NGC 4365) sug¬ 
gests that KDCs could often be integral parts of the com¬ 
plete bodies of galaxies. It is the balance of orbits that 
makes KDCs visible. More KDC galaxies, spanning a range 
of masses, should be observed and analysed in a similar way 
to determine if NGC 4365 and NGC 5813 are typical or spe¬ 
cial objects, and if the KDCs are decoupled or integrated 
structures. The spatial sampling, field-of-view, the spectral 
resolution and its efficiency make MUSE the right instru¬ 
ment to achieve this. 
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APPENDIX A: COMPARISON OF THE 
SCHWARZSCHILD DYNAMICAL MODEL AND 
THE MUSE KINEMATICS 

In Fig. lAll we show a comparison between the MUSE kine¬ 
matics and the Schwarzschild dynamical model from Sec¬ 
tion [5] The MUSE data were symmetrised using the point- 
(anti)symmetry as the model is by construction axisymmet- 
ric. The quality of the ht can be seen in the bottom row of 
maps which show the residuals between the model and the 
data, divided by the errors. The Schwarzschild model gen¬ 
erally reproduced well the kinematics, with an exception of 
the centre, which suggests the model likely requires a higher 
black hole mass. 


APPENDIX B: OTHER EMISSION-LINES 
DETECTED WITHIN THE MUSE 
WAVELENGTH RANGE 

As described in Section (6] we fitted a set of emission¬ 
lines typically found in the MUSE wavelength range: 
H/3, [Oiii]AA4959,5007, [Ni]AA5198, [Ni]AA5200, 

[Oi]AA6300,6364, Ha, [Nii]AA6547,6583 and [Sii]AA6716,673. 
Of these only All lines expect [Nl] is not detected, and [On] 
is detected on a much smaller area compared to the other 
lines. Ha and [Nil] are shown in the main text (Fig. [T] and 
in Fig. [BT] we show the distribution and kinematics of all 
other deteceted lines. 
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Figure Al. Comparison between the MUSE kinematics (top) and the Schwarzschild model kinematic predictions (middle) from Section[5] 
From left to right are shown: the mean velocity, the velocity dispersion and the higher-order Gauss-Hermite moments - Hq. The bottom 
row shows the residuals between the Schwarzschild model and the MUSE kinematics, divided by the observational errors. 
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Figure Bl. From top to bottom: Emission-lines H/3, [OIII|A4959,5007, [OI]A6300, [SII]A6716 and [SII]A6731. From left to right: A/rN 
map, logarithm of emission-line flux (units of of 10”^^ erg/s/cm^), the mean velocity and the velocity dispersion of a single Gaussian. 
Extraction of these emission-lines is explained in Section (6] but note that the velocity and velocity dispersion of H/3 was fixed to that 
of the Hq and of the forbidden lines to [NII|. This means that the lower four rows have the same kinematics, but different fluxes and 
distributions. Contours are the isophotes from the MUSE white light image. 
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